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Abstract 

We consider the spherically symmetric, asymptotically flat, non- 
vacuum Einstein equations, using as matter model a coUisionless gas 
as described by the Vlasov equation. We find explicit conditions on 
the initial data which guarantee the formation of a trapped surface in 
the evolution which in particular implies that weak cosmic censorship 
holds for these data. We also analyze the evolution of solutions after 
a trapped surface has formed and we show that the event horizon is 
future complete. Furthermore we find that the apparent horizon and 
the event horizon do not coincide. This behavior is analogous to what 
is found in certain Vaidya spacetimes. The analysis is carried out in 
Eddington-Finkelstein coordinates. 



1 Introduction 

A major open problem in the mathematical analysis of General Relativity are 
the geometric properties of spacetimes in the case when singularities develop 
out of regular initial data. Of particular interest is the validity of the cosmic 
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censorship conjecture which in non-technical terms states that genericahy 
every spacetime singularity which evolves out of regular, asymptotically flat 
data is covered by an event horizon so that it cannot be observed from infin- 
ity. For a precise form (or precise forms) of a statement of cosmic censorship 
we refer to [29] or [TT] . An important insight from the mathematical analysis 
of this problem is that the answer depends on the model used to describe 
the matter which collapses. For dust which by definition is a perfect fluid 
with pressure zero Christodoulou [7] showed that singularities which violate 
cosmic censorship — so-called naked singularities — arise for an open subset of 
spherically symmetric, asymptotically flat data. Christodoulou also carried 
out an extensive analysis of this problem with a massless scalar fleld as mat- 
ter model [8l[9l[T0]. Again, there do exist spherically symmetric data which 
give rise to naked singularities, but this behavior is shown to be unstable, 
and for generic data the spacetime has a complete future null infinity which 
is a more formal statement of the validity of the weak cosmic censorship 
conjecture, cf. [TT] . 

In [21 m [5] gravitational collapse was studied using as matter model a 
collisionless gas, i.e., the Vlasov equation. The Einstein- Vlasov system de- 
scribes an ensemble of particles which interact only through gravity and 
where collisions among the particles are neglected. It is a model for astro- 
physical systems like galaxies or globular clusters, although in astrophysics 
relativistic effects are most of the time neglected and the non-relativistic 
limit [22\ |26] of the system, the Vlasov-Poisson system, is used, cf. [6] and 
the references therein. The main result for the Einstein- Vlasov system in 
the cited investigations is that a class of spherically symmetric, asymptoti- 
cally fiat initial data is identified which lead to gravitational collapse and to 
a spacetime with a complete future null infinity so that cosmic censorship 
holds in the sense of [TT]. The analysis is carried out in Schwarzschild co- 
ordinates. It is shown that a piece of the maximal development of the data 
is covered by these coordinates which is sufficiently large to deduce the de- 
sired conclusions, but the analysis is hampered by the well known fact that 
Schwarzschild coordinates cannot cover regions of spacetime which contain 
trapped surfaces. Closely related to this defect is the fact that these coor- 
dinates do not seem to completely cover the radially outgoing null geodesic 
which generates the event horizon; at least it has not been possible to show 
that they do. 

It is hence natural to attempt an analogous analysis of the Einstein- 
Vlasov system using coordinates which do not break down at trapped sur- 
faces. In the present paper we propose to do so using Eddington-Finkelstein 
coordinates. For the purpose of illustration we recall that the Schwarzschild 
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metric, when written in Schwarzschild coordinates, takes the form 



Here M > is a constant, the time coordinate t G M coincides with the 
proper time of an observer who is at rest at spatial infinity, the area radius 
r > labels the surfaces of symmetry, i.e., the orbits of the group S0(3) 
which have surface area Airr'^ as measured in this metric, and 6 S [0,7r], 
G [0, 27r] parametrize these symmetry orbits. The metric is non-singular 
with signature ( — h H — h ) in the region r > 2M, but the coordinates break 
down at r = 2M. In Eddington-Finkelstein coordinates t is replaced by 

y ■-t + r + 2Mln(r - 2M). 

The Schwarzschild metric takes the form 

ds^ = ~ ~ ~~~) + '^^^ '^^ + ^^(^^^ + ^^^^ ^'^'^^) 

which is completely regular for r > 0, and r = is revealed as the true 
spacetime singularity where for example 

R"'^^^Raf3yS ~ r-^ as r ^ 0; 

^afd'y Riemann curvature tensor. The lines v = const correspond 

to radially ingoing null geodesies, and v is referred to as an advanced null 
coordinate or advanced time [Hi 5.5]. The line r = 2M is also null and 
represents the event horizon, while the surfaces of constant v and r are 
trapped if r < 2M in the sense that with increasing v also the radially 
outgoing null geodesies dr/dv = 1 — 2M/r at r < 2M move towards the 
center. We want to use such coordinates to study dynamic spacetimes so 
we generalize the above form of the Schwarzschild metric as follows: 

ds^ = -a{v, r) b^{v, r) dv"^ + 2b{v, r) dv dr + [dO'^ + sin^ 6 dip"^) . 

While b is required to be strictly positive, a can change sign, and as long 
as 6 > this metric is non-degenerate with signature ( — h -|- -|- ). Asymp- 
totic flatness means that the metric quantities a and b satisfy the boundary 
conditions 

lim a{v,r) = lim b{v,r) = 1. (1-1) 
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For a metric of this form the non-trivial components of the Einstein equa- 
tions 

Gap = 8TTTaf3 



are found to be 



(^rd^a + rabdra + a% - ab^ = SttTqo, (1.2) 
^ (rdra + a - l) = SttTqi, (1.3) 



\drb = 8TTTii, (1.4) 



rb 



(2drvb + b^drrd + 2abdrrb ^-r^ h Sbdrttdrb 

2o^ \ 



2b'^dra 2abdrb\ ^ _ 
+ — + = 87rr22; (1.5) 



the also non-trivial 33 component is a multiple of the 22 component. 

In a coUisionless gas the world lines of the particles are timelike geodesies. 
The ensemble can be described by a non-negative number density function 
/ on the tangent bundle TM or equivalently on the cotangent bundle TM* 
of the spacetime. The latter choice turns out to be advantageous for our 
analysis. We denote by {p^ , , , p^) the canonical momenta corresponding 
to the coordinates {x^ , , , x^) = {v,r,9,ip) so that {v,r,9,ip,po,pi,P2,P3) 
coordinatize the cotangent bundle TM*. On TM* the geodesic equations 
take the form 

dx"" a aB dpa 1 dg'^"/ 

where g'^f^ is the inverse of the Lorentz metric ga/3- Due to spherical sym- 
metry the angular momentum 

is conserved along geodesies, 

ar 

We shall as usual assume that all the particles in the ensemble have the same 
rest mass which we normalize to unity. Then their density / is supported 
on the mass shell defined by 

2 L 

-1 = g"^PaPf3 = rPoPi + a{pif + ^. 
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This implies that pi ^ always, and since we want dv/dr = = pi/b > 0, 
i.e., all particles move forward in advanced time, we require pi > and can 
express po as 

Po = -^(api + ^^j. (1.6) 

Due to spherical symmetry the particle density / is a function of the variables 
{v,r,pi, L). Rewriting the relevant components of the geodesic equations 
using V as the parameter we find the characteristic system (cf. ()2.7p . (|2.8p 
below) of the first order conservation law, i.e., the Vlasov equation, satisfied 
by/: 

+ ^ - ^r{a b)p, - drb ' 9p J = 0. (1.7) 

2 \r^px pi ) 

In order to close the system we have to define the energy momentum tensor 
in terms of / and the metric. In general, 

rj. I 1-1/2 /" dpodpidp2dp3 

Tap = \g\ ' I PaPfiJ , 

J m 

where \g\ denotes the modulus of the determinant of the metric and m the 
rest mass of the particle with coordinates (x^jPf^). In the above coordinates 
and using the restriction to the mass shell, 

/"OO ;"00 / \2 

Too{v,r) = - / S£0Lf{v,r,puL)dLdpu (1-8) 
Jo Jo Pi 



CO fOO 



Toi{v,r) = ^l I pof{v,r,pi,L)dLdpi, (1.9) 

Jo Jo 

Tii{v,r) = ^ / pif(v,r,pi,L)dLdpi, (1.10) 

f Jo Jo 

f'OO f'OO T 

T22{v,r) = -^ / -f{v,r,pi,L)dLdpi, (1.11) 
Jo Jo Pi 

where po has to be expressed via (jl.6|) . 

As our main result for the spherically symmetric and asymptotically flat 
Einstein- Vlasov system ()l.ip - ()l.lip we prove that data given at = which 
satisfy certain explicit conditions and do not contain a trapped surface, i.e., 
a(0, •) > 0, launch solutions where a{v,r) < for some t> > 0,r > 0, i.e., 
a trapped surface forms. A consequence of this result is that weak cosmic 
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censorship holds for our initial data in view of the results \12\ 113] . The 
solutions we construct have the additional property that for v sufficiently 
large, all the matter is strictly within {r < 2M} and the generator of the 
event horizon, which coincides with r = 2M for v large, is future complete. 
The structure of the data at t; = is essentially that a static, v independent 
state is surrounded by a shell of matter, and the particles in this shell move 
towards the center in a specified way. We also study solutions which at 
V = contain a black hole of mass < mo < M where M denotes the 
ADM mass for the complete spacetime. The data are such that the initial 
black hole is irradiated by Vlasov matter so that the apparent horizon is 
at r = 2mo for small advanced time v and then grows from r = 2mo to 
r = 2M in finite advanced time. The event horizon coincides with r = 2M 
for large v but lies strictly between r = 2mo and r = 2M for small v so 
in particular, apparent horizon and event horizon do not coincide. This 
behavior is analogous to what is found in certain Vaidya spacetimes where 
matter, however, is described by the ad-hoc model of null dust, cf. [171 5.1.8]. 

Our main result described above should be related to a previous study by 
Rendall [28] where he shows that there exist initial data for the spherically 
symmetric Einstein- Vlasov system such that trapped surfaces form in the 
evolution. However, the proof in [28] rests on a continuity argument and 
it is not possible to tell whether a given initial data set will evolve into a 
spacetime containing a trapped surface which is in contrast to the explicit 
conditions derived in the present paper. 

The paper proceeds as follows. In the next section we establish a local 
existence and uniqueness result for the initial value problem to the system 
(|l.ip - (jl.lip . together with a criterion which allows the solution to be ex- 
tended as long as certain quantities are controlled. In addition we collect 
some general information on the solutions to be used in the further analysis. 
In Section [3] we show that for particles in the outer shell pi is bounded both 
from above and away from zero as long as the shell stays away from the 
center. This result is used in Section H] to prove the formation of a trapped 
surface out of regular data. An essential step here is to show that if the 
particles in the outer shell move towards the center initially in a specified 
way they continue to do so for v > 0. In the last section we study the so- 
lutions after a trapped surface has formed and also obtain the Vaidya type 
spacetimes mentioned above. 

We conclude this introduction with some further references to the lit- 
erature. Background on the Einstein- Vlasov system and relativistic kinetic 
theory can be found in [TJ [27] . A distinguishing feature of the Vlasov matter 
model is that in the Newtonian case, i.e., for the Vlasov-Poisson system, no 
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gravitational collapse occurs and solutions exist globally in time, cf. |15tll6j. 
For more background on the Vlasov-Poisson system we refer to |20j . 



2 Local existence, conservation laws, and a-priori- 
bounds 

In this section we collect a number of observations, estimates, and a-priori- 
bounds which will play a role in what follows and which also lead to a local 
existence result and a continuation criterion. 

To begin with, let us assume that a non-negative, compactly supported 
number density f{v) £ C^(]0, oopx [0, oo[) is given at some instant v > 0. 
We want to show that at that instant the metric is then determined explicitly 
in terms of f{v). Firstly, the field equation (II. 4p can be integrated to yield 

6(f , r) = exp ^— 47r ^ 7]Tii{v,r]) dr]^ . (2-1) 

Since the formula (jl.lOp for Tn does not contain a metric coefficient, ()2.ip 
defines b{v) in terms of /(f), and 

< b{v,r) < 1, r > 0, (2.2) 

in particular b is positive as required. The boundary condition (II. ip holds 
as well, indeed b(v,r) = 1 when r lies to the right of supp/(t;). In order to 
express a we observe that 

Toi = -^{aTn + S) (2.3) 

where 

^ (2.4) 

Jo Jo Pi 
The field equation (jl.3p can now be rewritten in the form 

dr{ra -r)+ A-nrTu {ra - r) = -A-nr'^{Tii + S). (2.5) 

Since (ra — r)y^Q = this can be integrated to give 



a(f , r) 



1 - ^ ^ 47177^ (Til + S) exp ^- j A-KaTiida^ d-q; (2.6) 



clearly a{v,oo) = 1. Hence given f(v) at some instant v the metric is 
explicitly determined. Notice that when in Section [5] we consider the system 
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with a black hole at the center we replace the above boundary condition for 
ra — r at r = by one at r = oo. 

If the metric coefficients o and b are given and sufficiently smooth on 
some interval [0,V[ we denote by {R, Pi){s,v,r,pi, L) the solution of the 
characteristic system 

dr _b / l + L/r^\ 
^^"2^" {Pir J' 

dpi I (2hL i + L/r'^ \ 

— = - dr{ab)pi -drb 

ds I \r-^pi pi J 

of the Vlasov equation (jl.7p with (R, Pi){v,v,r,pi, L) = {r,pi). Then 

f{v,r,pi,L) = f{{R,Pi){0,v,r,pi,L),L) (2.9) 
is the solution of the Vlasov equation satisfying the initial condition /|„=o = 

o o 

/. If / is non-negative and compactly supported then these properties are 
inherited by f{v). 

The above observations allow for the following iterative scheme. If initial 

o o 

data / are given we define fo{v, r, pi , L) := /(r, pi,L). If /„ is given we define 
Tii^n and Sn by substituting /„ into the formulas (jl.lOp and (j2.4p . Next we 
define a„ and bn through ()2.1I) and (12. 6p . Finally, we obtain the next iterate 
fn+i via (j2.9p . using a„ and 6„ in (|2.7p . (|2.8p . This iterative scheme leads 
to the following local existence result. 

Theorem 2.1 Let f e Q(]0, oo[ X [0, cxd[) be compactly supported and non- 
negative. Then there exists a unique solution f G (^""^([O, y[x]0, oo[^x [0, oo[) 
o/ i/ie system with /|„=o = / where V > 0. Let V be chosen 

maximal. If 

sup I pi + — + - I {r,pi,L) G supp/(t'), 0<t'<l/l <oo 
{ Pi r J 

t/ien y = oo. 

Proof. In Schwarzschild coordinates the analogous proof has been carried 
out in detail in [19l[2T]. Here we only address some key issues. Firstly, by 
construction the iterative scheme indicated above will converge to a solution 
of the subsystem where only the field equations (jl.3p and (jl.4p hold. The 
remaining field equations can then be derived exploiting the following ob- 
servation. If Vq denotes the covariant derivative corresponding to the given 
metric, then the Vlasov equation implies that V^T"^ = 0. Also the Einstein 
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(2.7) 
(2.8) 



tensor satisfies the relation VqG"^ = 0. The relation VaE"' = satisfied 
by E"^^ := G"^ — SvrT"'^ then implies that the remaining field equations 
hold. 

Clearly, the supremum Q in the statement of the theorem together with 
()2.9|) controls Tn and S, hence a and 6, and via the corresponding field 
equations also their first order derivatives dra, drb, d^a; notice that 

- -d^a = Goo + abGoi = 8^(Too + aftToi). (2.10) 
r 

But in order to extend / as a solution it is necessary to control the 
derivatives of R and Pi with respect to their data r,pi,L which seems to 
require second order derivatives of a and b which are not controlled by the 
supremum Q. However, the quantities 

^ = drR, T] = drPl - Pl^drR 

satisfy a system of ordinary differential equations where second order deriva- 
tives of the metric coefficients appear only in the combination 

2drvb + b^drra + 2abdrrb 



which appears in the field equation ()1.5I) and is therefore controlled by Q as 
well; notice that d^b which appears in this combination cannot be controlled 
by itself in terms of Q. A similar argument also helps to obtain the bounds 
on the iterates required for their convergence. The differential geometric 
background of this maneuver is that the evolution of the derivatives of char- 
acteristics with respect to their data is governed by the geodesic deviation 
equation where derivatives of the metric enter only through the Riemann 
curvature tensor, and due to the symmetry the latter contains second order 
derivatives of the metric coefficients in the same combination as they appear 
in the Einstein tensor and hence in the field equations. □ 

It seems necessary that the support of the matter for the solution and 
hence for the initial data is bounded away from the origin. At first glance this 
is due to using polar coordinates (r, 6, 99), and one may hope to cure this by 
passing to the induced Cartesian coordinates r(sin0cos(/9, sin^sinc^jCos^), 
as was done in [19^ I21j . However, in the case of Eddington-Finkelstein 
coordinates this does not seem to help; the metric, when written in Cartesian 
coordinates, is not at the origin, indeed, not even continuous. Since this 
issue does not play a role in out analysis we do not pursue it further. 
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In what follows below, two situations will be considered. In the first we 
place a static, f-independent solution in the center and surround it by a 
shell of non-static Vlasov matter. Since we will only need to analyze this 
case as long as the outer non-static matter does not reach the central static 
part, the considerations above are sufficient to deal with it. In the second 
case, studied in Section [5l we consider the system only for r > where 
a(0, ro) < 0, so again a neighborhood of the center is avoided. 

Coming back to the first case we need to see that the system indeed has 
static, ?;-independent solutions. The existence of static, spherically symmet- 
ric and compactly supported solutions has been established in |18i [23l [M] us- 
ing Schwarzschild coordinates. Given such a static solution in Schwarzschild 
coordinates its metric takes the form 

ds^ = -e't^(-)dt^ + e^^'^'-Ur' + r^dO^ + sin^ 6 d^'). 

Changing to the variable v = v{t, r) = t + e^^'^'^''^^'^'^ dr] brings this metric 
into the Eddington-Finkelstein form where 

a = e-^\ b = e^+^ 

so that these transformed metric coefficients are independent ofv as desired, 
and the same is then true for the particle density / of the steady state. 

In what follows fs denotes a fixed spherically symmetric steady state 
of the Einstein- Vlasov system with spatial support in [0, rg]. We consider 
initial data of the form 

o o o o 

/:= /s + /out where /out G C'^ (]r-o,oo[x]0,oo[x[0,oo[), /out > 0. (2.11) 

For a solution launched by such data we introduce the notation 

fout{v,r,p^,L) := /out((^,^'i)(0,i^,r,pi,L),L), (2.12) 

i-e., /out is the density of the particles which initially do not belong to the 
steady state. Due to spherical symmetry there are no gravitational waves, 
and the only way that /out can influence the matter in the region {r < tq} is 
if outer particles actually reach this region. As long as this does not happen 
the matter in the central region stays in its equilibrium configuration. Hence 
we obtain the following corollary to Theorem 12. 1[ 

o _ 

Corollary 2.2 Initial data f as specified in Ii2.11\) launch a unique solution 
f £ Cn[0,F[x]0,oopx[0,oo[) of the system flJl i^/ fOT]) with /|„=o = / 
where V > 0. Let V be chosen maximal. If 

suplpi + — I (r,pi,L) G supp/out(^^), < < fI < oo 
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and 

inf {r I (r,pi,L) G supp /out(t^), < v < V} > vq 

then V = oo. 

We now collect some properties of the local solutions which are used in what 
follows. First we look at the ADM mass and resulting a-priori bounds. 

Proposition 2.3 Let 

{ \ (^ ( \\ ■ ( \ ^ 2m{v,r) 

m(v,r) :=-(! — a[v,r)) , i.e., a(v,r) = 1 



2 r 
Then the quasi-local ADM mass m is given by 

m{v,r) = 2iT j r/^(Tii + 5) exp ^— y ATTaTudcj^ di] 

= 27r ^\\T^^ + S){v,v)^-^dr^. (2.13) 
Jo b{v,r) 

For r sufficiently large, 

m{v, r) = M = lim m(u, r) 

1 — >oo 

which is a conserved quantity, the ADM mass. Moreover, 

m(t;,r) =27r/ r]^{aTii + S){v,rj) dr], (2.14) 



and 

a{v, r) < 1, 
b{v , r) m{v , r) = 2Tr / rf {Tn + S)b drj 



poo 

<2Tr rf{Tii + S)bdr] = M. (2.15) 

JO 

Proof. The formula (j2.13p for m{v,r) results from comparing its relation 
to a with (g^l) and ([21]). Recalling (I21IID . 

d^m = --d^a = 47rr2 ( Itoo + aToi ) (2.16) 



in particular, m(-,r) is constant for r large enough and hence to the right 
of the support of the matter. If for the moment we denote the right hand 
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side of (j2.14p by m then clearly rh{v,0) = = m{v,0) , and the differential 
equation (j2.5p for ra — r is equivalent to the fact that drfh = drm. By (j2.13p . 
m > and hence a < 1, and (j2.13p together with the conservation of M and 
the fact that b{v,oo) = 1 implies (|2.15p . □ 

To conclude this section we collect some information on sign changes in a. 

Proposition 2.4 //a(fo,To) < for some ^ 0,ro > then a(u,ro) < 
for all V > for which the solution exists. If a(v,r) < then all timelike or 
null geodesies at the spacetime point (v, r, 9, (p) move towards strictly smaller 
values of r. 

Proof. First we note that by ()2.16p . 
dva = -Svrr ^^^bo + aToij 

9_2l poo poo 1 . . 

Hence dva{v,r) = if f{v,r,-) = 0. If f{v,r,-) > then dya{v,r) < 0, 
provided a{v,r) is sufficiently close to so that the term in parentheses is 
positive; notice that f{v,r,-) has compact support. Hence if a{vo,ro) < 
at some vq , rg this sign must be preserved for all v > vq. The remaining 
assertion follows from the geodesic equations. □ 



3 A lower and an upper bound on pi 

In this section we consider solutions launched by data as specified in (j2.1ip . 
i.e., with a steady state at the center. The aim is to prove that pi is bounded 
from above and below on the support of /out where we recall (j2.12p . In 
connection with the continuation criterion in Corollary 12.21 this means that 
a solution can only blow up if particles from the exterior mass shell reach the 
central region {r < tq} where the steady state part is supported. The result 
may be compared to |25j where it is shown that in Schwarzschild coordinates 
solutions which blow up at all must do so at the center first. 

Theorem 3.1 Let V > be such that r > tq for (r,pi,L) € supp/out('w) 
and < V < V . Then there are constants Ci, C2 > such that Ci < pi < C2 
for (r,pi,L) G supp/(u) and < v < V . 

In the proof of this theorem we need the following a-priori-bound. 
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Lemma 3.2 Let y > and rg > 6e such that f{v, ro, •, •) = for < v < 
V. Then 

roo r-oo 

/ r'^Tn{v,r)dr < r^Tu{0,r)dr^ v,0<v<V. 

Jro Jro 27rro 

Proof. Using the Vlasov equation, integration by parts, and the field equa- 
tions we find that 

— / r ±11 ar = —- I brlii[\ — a)ar + ii Iff — jaLapidr. 

dv Jro 2 Jro Jo Jo Pi 

Now a < 1 and Tu > so the first term can be dropped. As to the second, 
we observe that 

bL l + L/r^l 
r^pi pi r 

and recalhng the expression (|2.4p for S we obtain the estimate 

^ /'oo POO ^ roo 

— / r^Tii dr < rbSdr < — / r'^bSdr. 

av Jro Jra ^0 Jro 



dv 

By (EIS]), 

poo poo 

M = 2-K I r'^{Tii + S) b{v, r) dr > 2-k / r'^S{v, r) b{v, r) dr, 

Jo Jro 
^^^^^ ^ oo 

— / r^Tii < 



dv Jro 27rro 
and the claim of the lemma follows. □ 
Proof of Theorem 13.11 Let [0,V] 3 s ^ {r{s),pi{s), L) denote a charac- 
teristic in supp/out- Using (j2.7p we can rewrite (|2.8p in the form 

In order to obtain a lower bound for pi we observe that 

+ -Grab + aOrb ; — — — 



ds Pi r'^{pi)^ \2 b ds J pi 

<r (^^ ^ K drhdr\ 1 

< -drab + adrb — — — . 

\2 b ds J pi 
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Consider an arbitrary instant of advanced time < vq < V. Applying 
Gronwall's lemma we find that we need to estimate the integral 



VQ 



^drab + adrbj {s,r{s)) - (^j {s,r{s))^{s) 



ds 



which is a curve integral along the curve 

7 = {(s,r(s)) I < s < vq}. 

In order to estimate this integral we apply Green's formula in the plane. For 
this we define 



Ci = {{vo,r) \r{vo)<r<R}, 
C2 = {{v,R) \ 0<v<vo}, 
C3 = {(0,r) I r(0) <r <R}. 

Here i? > is large, and we will let R ^ 00. We orient the closed curve 
r = 7 + Ci + C2 + C3 clockwise. Now 



VO 



-drab + adrb — — 

2 b as 



ds 



-drob + adrb \ dv — (—j^ 1 dr 



Ci JC2 JC3 

We denote by O the domain enclosed by T and apply Green's formula in the 
plane to find that 



1 



dr dv 



,drb drttb 



dr dv: 



recall the field equation (jl.Sp for the 22-component G22 of the Einstein 
tensor. If we now use the field equations and the definitions of the energy 
momentum tensor it turns out that 



22 



adrb drab 4tt^ f°° f'^l + 2L/r' 



r 

Hence by ^U^ . 



< 



r Jo Jo 
2mb 



Pi 



fdLdpi 



b 2m 



r 



n r' 



2 dr du < —vq- 
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Next we use the field equation (jl.4p and Lemma [3?2] to get the estimate 



/ •••=/' '-r^dr<4. rrTnivo,r)dr 
Jci Jr(vo) b{vo,r) 



< 



oo 



< — 47r / r''Tn{vo,r)dr (3.2) 



2r 

/ 



— Uvr / r2rii 0,r dr + t;o • 

The C2-contribution vanishes in the limit R ^ oo: 

r r° fi \ M 

-J ■■■ = y {^■^drab + adrbj{v,R)dv<^vo, 

since for r sufficiently large and in particular outside the support of /, 6 = 1 
and dra = 2M/r^. Finally, 



-/ ■■■=-/ 

Jc3 Jrio) KO,r) 

Altogether this implies, using Gronwall's lemma, the estimate 

1 1 /"°° 2^ . X , 3M \ 

— — <— — exp — / r Tu{0,r)dr + —V ] 
pi{v) pi{0) \ro Jro J 



(3.3) 



In order to estimate pi from above we start from (j3.ip . The term (bL) / (r'^pi) 
is bounded by the previous estimates so after using Gronwall's lemma we 
have to estimate the term 

-drab + adrbj dv — 1 dr 



+ / + 



from above, where the curves are defined and oriented as before. The Ci- 
contribution is now negative and can be dropped. Since for r sufficiently 
large and outside the support of /, 5 = 1 and dra >^ the C2-contribution 
can be dropped as well. The Ca-contribution is determined by the initial 
data, so it remains to estimate the integral over the closed curve F which we 
again turn into an integral over Q, using Green's formula. Due to the change 
in sign we are left with estimating the integral 

— T^b [ [ — — — fdLdpidrdv. 

n Jo Jo Pi 
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Using the already established bounds this amounts to estimating the inte- 
gral Jq°° / dL dpi dr which represents the number of particles in the 
domain {r > tq} and is conserved for v G [0, V]. The proof of Theorem 13.11 
is complete. □ 
Together with the continuation criterion from the local existence result 
Corollary 12.21 we obtain the following corollary. 

Corollary 3.3 There exists V > such that the solution exists on the in- 
terval [0, V[ and r > tq for all particles in the outer matter. IfV is chosen 
maximal, then V = oo if 

inf{r I (r,pi,L) G supp /out(t^), < v <V} > ro, 

i.e., the solution can be extended as long as the outer matter stays outside 
{r < ro}. 



4 The formation of a trapped surface 

In this section we want to specify conditions on the data such that a trapped 
surface evolves, but is not already present in the data. The data are again 
of the form (j2.1ip . i.e., they have a steady state fs at the center whose mass 
we denote by mg. We fix some notation: 

supp/out C [Ro,Ri] X [p-,p+] X [0,L+] (4.1) 

where 

< ro < Ro < Ri, < p_ < p+, L+ > 0, 

o o 

and M > mo is the mass of / := + /out- If ^ 2M then the data do 
not contain a trapped surface in the sense that a(0, •) > 0. We define 

P := -maxjp^ | {r,pi,L) E supp/out} • (4.2) 

Since 

1 1 1 11 + L/r2 1 11 ^ ^ 

P = tPo + api = -api - < -pi - -— (4.3) 

2 2 pi 2 2 pi 

we have P > 0, i.e., all particles move inward initially, if for example p+ < 1. 
The following theorem is the main result of the present paper. 



Theorem 4.1 Let data fg + /out be given such that 



Ro>2M> ro, L+ := 12m^, 2P > 1 + (4.4) 



r, 



2 ' 
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and such that there exists V > with the property that 

— — — ^exp )v>Ri-2M, (4.5) 

o 

Then the solution launched by fg + f forms a trapped surface at some ad- 
vanced time V < V . 

Remark. Note that this result impUes that weak cosmic censorship holds 
for these data in view of the results |12lll3j . In Section [5] we will see that for 
data as considered above an event horizon evolves which is future complete, 
cf. Theorem 15.21 and Corollarv l5.3i 

Remark. It is easy to see that data which satisfy the conditions above 
do exist: First we arbitrarily fix the central steady state fg and hence also 
ro and mo. Next we fix Rq,M,p^ in such a way that ()4.4|) holds; notice 
that P becomes large if we chose p+ > small, cf. (j4.3p . Then we choose 
< p_ < If we replace the inequality in (14. 6p by an equality this 
uniquely determines V which we then choose slightly smaller to preserve the 
inequality. Substituting this V into (j4.5p we see that that relation is satisfied 
provided Ri is sufficiently close to 2M, and all the support parameters and 
the mass of /out are fixed. Now we fix any non-negative, non-vanishing 
function g which satisfies the support conditions, consider data Ag with 
some amplitude ^ > and denote the induced ADM mass and metric 
coefficient by and ua respectively. From ()2.13p we see that Ma depends 
continuously on A with Mq = ttiq. If a a remains positive as A increases then 
()2.14p shows that Ma becomes as large as we wish with increasing A. So 
assume that aA*{r*) = for some r* > Rq and some value of amplitude A* , 
and OA is positive for all A < A* . Then 2mA* (r*) = r* > Rq > 2M and aA* 
is still non-negative, hence Ma* > M, and Ma = M for some smaller value 

o 

of A. Hence /out satisfying all support conditions and having the proper 
mass exist. Since Rq > 2M, a(0, •) > for any such /out- 

A major step in the proof of Theorem 14.11 is to show that, for suitable 
data, the particles move inward as long as they stay outside r = tq. Since 
we will recycle this result in Section [U we formulate it in a more general 
setting than needed here. 

Lemma 4.2 Let V be the length of the maximal existence interval of a 
solution f as in Corollary \3. 31 and assume that 

L+ < I2ml and P > 0. 
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Then 



max{p^ I {r,pi,L) £ supp/out('y), <v <V} < -P < 0. 

Proof. By continuity and the assumption on the data there is some interval 
[0, V*[C [0, V[ such that < for all (r,pi, L) G supp/out('v) with <v < 
V*. We choose V* maximal and have to show that V* = V. To this end 
we consider a characteristic in supp /out and parametrize it by proper time. 
After some computation we find that 



P 



d 



1 



-po + api 



drb 2 , dva 2 , "-drU 
--^{Po) +^iPi) +^ 



I n2 , ^r-a 



L 

PoPi + a— 



47rr 2 f 87rr 1 \ 

--^TniPo) + [j^^oi + ^(1 - a) j popi 

+ (— 5r^oo + ^(l-a)J {pif + a^ 



- a)PoPi + 7^(1 - a){pif + 
ro Zr 



47r^ 



oo roc 



JO 



Pi(Po)^ - 2popopi + 



{Po)' 

Pi 



-{Pl? 



f di dpi 



where pq is defined as in (jl.6p but in terms of the integration variables pi 
and L. Now 

PoPi ' 



PiPo 



iPi 



>0, 



so after substituting the definition (|1.6p for 



P^<-^(i-«) 1 + ^:2 



(1 a) ^ + ^ 



If a(u, r) < 0, then 



(4.7) 



1 ^ ( ^ 



L L 

^ + ^ < 0. 



If o(f,r) > 0, then we observe the relation a(f,r) = 1 — 2m[v,r) /r and 
investigate the behavior of m(f , r) with respect to u . By Proposition 12.41 
a(-,r) > on [0,^]. By (f2l6]) and ([THD, (fTa]) . 

■l(po)2 



= 47r 



JO 



6 



apo fdLdpi, 
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and 

l (Po)^ , Po 1 
h apQ = —p . 

b pi pi 

By assumption, p-*^ < on [0,v] C [0, V*[ for all particles, also pi > for all 
particles, and 

6/ 1 + 



2 V Pi J 

because of the sign of a(-,r) on [0,v]. Hence dvm{-,r) > on [0, f], and 
m{v,r) > m{0,r) > niQ. Hence (|4.7p implies that 

1 1 2m 

P" < 

2,r r 

< 1 (Lr - 3Lmo - r'ruo) = -J ^ - (v^or " ^ ) " 3Lmo 





\ 2m L 




)-—-,+ 




2 N 1 


'0 — 





< . ^ fL- 12m^) < 0. 




4r^mo 

Together with the case a{v,r) < this shows that p^ < for any char- 
acteristic in supp/out(v) with v £]0,V*[, hence V* = V, and the proof is 
complete. □ 

Proof of Theorem 14.11 Assume that a{v, r) > for all r > and < 

V < V and as long as the solution exists. We will show that the solution 
must then exist on the interval [0, V] and m{V, r) = M for some r < 2M so 
that a{y, r) < which is the desired contradiction. The basic idea is that 
(|4.6|) guarantees that the continuation criterion from Corollary 13.31 applies 
to the interval [0, V] so that the solution exists there, while (|4.5p guarantees 
that all the matter arrives inside {r < 2M} within that interval of advanced 
time. 

We first need to improve the a-priori-bound from Lemma 13.21 under the 
assumption that a > 0. By Lemma 14.21 and ()4.3p . 

1 ll + L/r^ . 

- o =P <-P 

2 2 pi 

on supp/out(w) for all v > 0; P > by (|4.4p . The fact that a > and the 
assumptions on the support of the initial data imply that 

Pi < < 1 < — < — + api (4.8) 

2-P pi pi 
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on supp/out(w). Since a > 0, (jM]), (EHI), and (|2J5]) imply that 



r^Tii{v,r)dr < I r"^ {S + aTu) {v,r) dr < —M. (4.9) 



ro 



ro 



2tt 



Along any characteristic in supp/out and for < v < V , since a > 0, 



6 / 1 + L/r^ 
2 



1 



a < - 1 + ^ 



1 



(Pi)^' 



Now we use the estimate (I3.3p . but we control the Tn integral in (|3.2p by 
Col) instead of Lemma |3.2[ Hence 



dr 
dv 



1 / L 



1 /2M^^ 4M 

exp V H 

^-0 



The assumption (14.61) now implies that r > vq for all particles in supp/out(v) 
and < V < V . Hence by Corollary 13.31 the solution exists on the interval 
[0,V]. 

We now show that at v = V all the matter must be strictly within 
r < 2M so that a{V, r) < for some r < 2M, which is a contradiction. 
First we note that by (14. 9p . 



r) = exp ^— 47r J r]Tii{v , t]) di]^ 

j ri^Tii{v,ri) dt]^ > exp 



> exp 



47r 



2M 



Hence along any characteristic in supp /out by (|4.8 



dr 
dv 



pi pi 

b > exp 



p 



,0 



Pi 



2M\ P 



ro J Pi 



> exp 



2M\ 2P2 



2 • 



Assumption ()4.5p now implies that by u = F all characteristics starting in 

o 

supp /out are strictly within r < 2M, and the proof is complete. □ 



5 After the formation of a trapped surface 

In this section we investigate the Einstein- Vlasov system (jl.ip - ljl.lip for 
data which are such that a(0, ro) < for some ro > 0, i.e., the data already 
contain a trapped surface. Besides being of interest in itself this problem is 
relevant in connection with the result of the previous section, where it was 
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shown that such data evolve out of data not containing a trapped surface, 
but where the argument stopped after a becomes negative somewhere; notice 
that Theorem 14.11 was proven by contradiction. If we for example want to 
show that eventually all matter ends up in the region {r < 2M} and that 
r = 2M is complete and the event horizon of the evolving black hole we have 
to be able to continue the analysis after a has become negative somewhere. 

First we check that the system is well posed on the domain {r > ro}. 
If / and therefore Tu and S are given for r > ro, we define b by (j2.ip as 
before. Next we recall (12. 5p . The desired boundary condition for ra — r at 
infinity is liuir^ooifa — r) = —2M, and hence 

2m{v,r) 

a{v,r) = l , (5.1) 



where 



m{v,r)=M-^j^ 4Trr]'^{Tu + S)exp(^- j ATvaTuda 



dr] (5.2) 



and M > denotes the ADM mass of the spacetime. 

If a(0,ro) < for some tq > then by Proposition 12.41 O'ivjro) < for 
all w > for which the solution exists. Since characteristics can only leave 
but never enter the region {r > tq} when followed forward in v, 

fiv,r,p^,L) = /((i?,Pi,L)(0,^;,r,pi,L)) 

defines the solution of the Vlasov equation ()1.7p on {r > ro} with initial 

o 

data /. 

We now specify the data which we consider in the present section: 
/ S C-'^([ro, oo[x]0, oo[x [0, oo[) is non-negative and compactly supported 
in [ro, oo[x]0, oo[x [0, oo[. Defining 

f°° 2 ° ° f n ° \ 

mont := 27r / rj {Tu + S) exp ( - / AnaTuda 1 drj 
we choose some constant 

M > mo,t + y 



which plays the role of the ADM mass in (|5.2p . Then 

2m(0, ro) 
a(0,ro) = 1 < 
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as desired. Notice that the system can be studied on {r > tq} regardless 
of what represents the mass in the central region {r < rg} as long as that 
mass is there initially. In particular, this has the advantage that we need 
not take care that the outer matter does not interfere with a central steady 
state as was necessary in Theorem 14. II First, we obtain the following global 
existence result. 

o 

Theorem 5.1 Initial data f as specified above launch a unique solution 
f £ C"'^([0, oo[x [ro, oo[x]0, oo[x [0, oo[) of the Einstein- Vlasov system il.l]) - 

Proof. As a first step one can establish a local existence result correspond- 
ing to Theorem 12. H using the arguments indicated there and in the well- 
posedness discussion above. If F > is the length of the maximal existence 
interval and 

sup <pi H I (r,pi,L) e supp/(u), < V < v \ <oo 

I Pi J 

then V = oo. But the required lower and upper bounds on pi along charac- 
teristics follow by the estimates in the proof of Theorem 13.11 □ 
Let us now define mo := M — mout > as the mass initially in the region 
]0,ro], and assume that 

supp/ C [ro,Ri] X [p-,p+] X [0,L+] 

where 

Ri > ro, < p_ < p+, L+ := 12m^, 
and we require that all particles move inward initially, more precisely. 



2F>l+p^. (5.3) 

where as before, 

P := — max | (r,pi, L) G supp /| . 
We obtain the following asymptotic behavior for large advanced time. 

o 

Theorem 5.2 For data f as specified above all particles in the region {r > 
ro} are moving inwards, p^ < —P. Moreover, there exist V* > and 
ro < R* < 2M such that for v >V* all the matter is in the region {r < R*}, 
and a{v,2M) = 0. The line v > V*, r = 2M is a radially outgoing null 
geodesic which is future complete and is the generator of the event horizon 
of the spacetime for v >V* . 
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Proof. First we observe that we can apply Lemma 14.21 also in the present 
situation on the region {r > tq}; notice that in the proof of that lemma no 
assumption on the sign of a was made. Hence < —P for all the particles 
in that region. This implies that 



dr 
dv 



Pi Pi 



for these particles. As long as r > 2M and hence a > it follows that 
h > exp(— 2M/2M) = 1/e; for this estimate we again rely on (j4.9p so we 
need the estimate corresponding to ()4.8p . but now only for r > 2M which 
accounts for the relaxed condition (j5.3p . By (j4.3p and since a > for 
r > 2M, 

= -2pi + api > 2P 



so that 



dr 
dv 



2p2 

> ^ > 



e(l + L/r2) - e(l + L+/(2M)2) 

as long as r > 2M. This shows that there exists > such that for v >Vi 
all particles are in the region {r < 2M}. If we choose some ro < R < 2M 
and V* := Vi + 1, then 

c := inf 1^^^ \ Vi<v<V*, R<r<2M, (r,pi,L) G supp/(?;)| > 0. 

This implies that within the interval [Vi, V*] all particles must have moved 
by a uniform distance to the left which proves the existence of R* . 
For v>V* andr> R*, 

a{v,r) = 1 - 2M/r, b{v,r) = 1. 

The line r = 2M, v > V* is null since a{v, 2M) = 0, and examining the 
geodesic equation for this null geodesic shows that it exists on an interval 
of affine parameter which is unbounded to the right. The proof is complete. 

□ 

If we go back to Theorem 14.11 we see that for a solution as considered in 
that theorem there exists v > such that fontiv) defines data which satisfy 
the assumptions of Theorem 15.21 In particular, there exists < ro < 2M 
such that a(f,ro) < 0, and we obtain the following corollary. 



Corollary 5.3 The solutions obtained in Theorem \4-l\ exhibit the same 
asymptotic properties as obtained in Theorem\5^ 
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To conclude this paper we exploit our estimates to construct a class of 
solutions where the initial data represent a black hole surrounded by a shell 
of Vlasov matter. These solutions will illustrate the fact that event horizons 
and apparent horizons do in general not coincide. 

Remark. The fact that in general the apparent horizon and the event 
horizon need not coincide is usually illustrated by sending a shell of so- 
called null dust, i.e., a pressure-less fluid of photons, into a black hole, 
and the corresponding spacetimes are known as Vaidya spacetimes. The 
theorem below shows that the corresponding behavior of the horizons can 
also be achieved with a less artificial matter model. 

o 

We consider non-negative and compactly supported data / G 
C^(]0, oopx [0, oo[) with mass 

j-oo / poo \ 

< mout = y rf{f II + S) exp [- J Anafu daj dr] < M 

and 

supp/ C [-Ro,^i] X [P-,P+] X [0,L+] 

where 

< 2M < i?o < ^1, 0<p- < p+, L+ = 12ml = 12(Af - mout)^- 

As before, we require that all particles are initially moving inward sufficiently 
fast in the sense that (j5.3p holds. For r* := 2mo < 2M we have a(0, r*) = 
while a(0, r*) > 0/ < for r > r* /r < r* . In other words we have a trapped 
region < r < r*, a black hole of mass mo, surrounded by a shell of Vlasov 
matter which is moving inwards, and we obtain the following result. 

Theorem 5.4 Data as specified above launch a unique solution f on the 
domain [0, oo[x]0, oo[. All the particles in supp/ are moving towards the 
center which they all reach within a finite interval [0, r*] of proper time. 
There exist V* > and < R* < 2M such that m{v, R*) = M for v > V* . 
Hence if 

R{v) := sup{r > | a{v,r) < 0} 

then R{v) = r* on some interval of advanced time [0, v*] with < v* <V*, 
and R{v) = 2M for v >V* . On the other hand, the generator of the event 
horizon is a radially outgoing null geodesic which coincides with R{v) = 2M 
for V > V* , but lies strictly to the right of R{v) = r* for v G [0, f*]. 
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Proof. If we take any < tq < r* we can apply Theorem 15.11 to obtain 
a solution on [0, oo[x]ro, oo[. If we decrease rg we get an extension of this 
solution, and since < tq < r* can be arbitrary the have the solution on 
the asserted domain [0, oo[x]0, oo[. 

Next we apply Lemma 14.21 which shows that all particles continue to 
move inwards with sup{p^ | {v,r,pi, L) £ supp/} < which proves the 
assertion on their behavior in proper time. That all particles end up strictly 
inside {r < 2M} within a finite interval [0, V*] of advanced time can be 
shown exactly as above, and the assertions on the event horizon and the 
apparent horizon follow. □ 

Concluding Remark. The analysis in the present paper leaves open the 
question of what happens at the center r = 0. Numerical evidence suggests 
that in the situation of Theorem 14.11 respectively Theorem 15.21 a spacetime 
singularity arises where the Kretschmann scalar R"'^^^ Raises blows up at 
r = 0. At the same time the region where a < extends all the way to r = 
so that no causal curve can enter {r > 0} out of the singularity and strong 
cosmic censorship holds. It remains to prove these assertions. 
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